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ABSTRACT 

Several groups have recently computed the gravitational radiation recoil produced by the merger of 
two spinning black holes. The results suggest that spin can be the dominant contributor to the kick, 
with reported recoil speeds of hundreds to even thousands of kilometers per second. The parameter 
space of spin kicks is large, however, and it is ultimately desirable to have a simple formula that gives 
the approximate magnitude of the kick given a mass ratio, spin magnitudes, and spin orientations. 
As a step toward this goal, we perform a systematic study of the recoil speeds from mergers of black 
holes with mass ratio q = toi/to2 = 2/3 and dimensionless spin parameters of oi/toi and a2/m2 equal 
to or 0.2, with directions aligned or anti- aligned with the orbital angular momentum. We also run 
an equal-mass ai/mi = —a2/m2 = 0.2 case, and find good agreement with previous results. We find 
that, for currently reported kicks from aligned or anti- aligned spins, a simple kick formula inspired by 
post-Newtonian analyses can reproduce the numerical results to better than ~10%. 
Subject headings: black hole physics - galaxies: nuclei - gravitational waves — relativity 



1. INTRODUCTION 

In the past two years, numerical relativity has under- 
gone a revolution such that now multiple codes are able 
to stably evolve the last few cycles of the inspiral, nierger , 
and ringdown of two black holes ("Briigmann et al .ll2004l: 
Prctorius 2005; CampancUi ct al. 2006a, b; Bake r et al.l 
2006a| bl). An important astrophysical output of 
such simulations is the net recoil due to asymmet- 
ric emission of gravitational radiation, because this 
has major implications for the growth of super- 
massive black holes in hierarchical merger scenarios 



(iMerritt et al.ll2004: Bovlan-Kolchin ct al. 2004: Haiman 
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holes and curre n t-day i ntermediate-mass b la ck hol es 
(iTaniguchi et all l2000t iMiller fc H amilton 20023,1?; 



Mouri fc Taniguchil l2002al l bl: iMiller fc Colbert 200 
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BekensteinI 

^ iFitchett fc Dctwcilcr 1984: 

Redmount fc ReesI '1989': Wiseman 1992; Favata et 
2004 : iBlanchet et al. ,^005; Damour fc Gopakuma: 
20061 ) have now been augmented with full numerical 
results for non-spinning black holes with different 
mass ratios (iHe rrmann ct al. 2007c; Bak er et al.ll2006"3 
IGonzalez et al.li2007al) , for black holes with equal masses 
and spins initially orthogonal to the orb ital plane 
(jHerrmann et al.ll2007at iKoppitz et alll2007t l. for black 
holes with eq ual masses and spins init i ally parallel to the 
orbital plane (|Gonzalez et all l2007bt ICampanelli et all 
l2007bV for black holes with equal masses and spins 
initially oriented at some angle between the orbital plane 
and the orbital angular mome ntum (|Herrmann et all 
I2007bt iTichv fc Marronettil [200l . and for black holes 
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with unequal masses and spins initially either parallel 
to the orbital plane or oriented at some angle between 
the orbital plane and t he orbital angular momentum 
(jCampanem et al.|[2007al ). 

The parameter space for kicks with spin is large, so 
for astrophysical purposes it is important to have a sim- 
ple parameterized formula for the kick that can be in- 
cluded in simulations of N-body dynamics or cluster and 
galaxy m ergers. For non-spinning black holes, the classic 
IFitchettI ([1983 ) formula V oc rj^ {5mlM) does a reasona ble 
job (although not perfect; see IGonzalez et all l2007al for 
numerical results), where for black hole masses mi and 
7712 > TOi we define M = mi + 7772, 5m — 7712 — 7711, and 
77 = 777irr72/M^. Initial explorations of equal-mass spin 
kicks also show evidence of sim plicity, with fits li near in 
spi n reproducing the res ults of IHerrmann et all (|2007af l 
andl Koppitz et all (|2007f ). This encourages us to explore 
a more general class of kicks. 

Spins that are aligned (prograde) with the orbital 
angular momentum may be of particular interest be- 
cause it has been argued that interaction with accre- 
tion disks will tend to a lign spins during "wet" mergers 
(jBogdanovic et al.ll2007| ). The only previous numerical 
studies for a black hole with an ali gned spin were car- 
ried g ut by IHerrmann et all (|2007af) and iKoppitz et all 
()2007D . which considered only cases where the spin of the 
second black hole was anti-aligned (retrograde) with the 
angular momentum and the masses of the black holes 
were equal. Here we compute the kick speeds from a 
q = 7771/7772 = 2/3 mass-ratio set of mergers, with spin 
parameters of or 0.2, and directions either ahgned (pro- 
grade) or anti-aligned (retrograde) with the orbital an- 
gular momentum. The symmetry of the configuration 
therefore guarantees that the kick direction is in the or- 
bital plane. We find a fo rmula that matches all o f our 
kick speeds, and tho se of IHerrmann et ah! (j2007af ) and 
IKoppitz etall (|2007f ). to within ~ 10%. If the in-plane 
kicks can be generalized straightforwardly to more gen- 
eral orientations, and added to kicks perpendicular to 
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the orbital plane, there is the prospect of simple astro- 
physical modeling of the gravitational rocket effect for 
arbitrary black hole mergers. In § 2 we describe our ini- 
tial data and methodology. In § 3 we present our results, 
and discuss the implications of these simulations. 

2. INITIAL DATA AND METHODOLOGY 

In the following, we use geometrized units where New- 
ton's gravitational constant G and the speed of light c 
are set to unity so that all relevant quantities can be rep- 
resented in terms of their mass-scaling. For example, 1 
Mq is equivalent to a distance of 1.4771 x lO^cm, and a 
time of 4.9272 x 10~^s. Accordingly we express distances 
in terms of M, the initial (ADM) mass of the system. 

We simulated inspiraling black-hole binaries of various 
mass ratios and spins, with the same initial coordinate 
separation in each case. In these cases our initial 
mass ratio approximated either 2/3 or unity. Our 
simulations wer e performed with our finite differencing 
code Hahndol (jimbiriba et al.l l2004l ). which solves a 
3-t-l formulation of Einstein's equations. Adaptive 
mesh refinement and most parallelizati on was handled 
by the software package Param esh (IMacNeice et aL 
[2000) . For initial data we used the iBrandt fc Briigmann 



(fr997i) Cauchy surface for black hole punctures, as 
computed by the second-order-accurat e, multigrid 
elHptic solver AMRMG (jBrown fc Lowdl2005l ). We evolved 
this data with the standard Baumgarte-Shapiro- 
Shibata-Nakamura (iNakamura. Oohara fc Koiimal 



1987'; 'Shibata fc Nakamura"1995'; 'Baumgarte & Shapirc 
1998; Imbiriba ct al. 2004) evolution equa ,tionSj modified 
only slightly with dissipation terms as inlHiibned (Il999l) 
and constraint-damping terms as in iDuez et al.l (|2004D . 
Our gauge conditions were chosen a ccording to the 
"movi ng puncture" method, as in Ivan Meter et al.l 
(|2006f ). Time-integration was performed with a fourth- 
order Runge-Kutta algorithm, and spatial differencing 
with fourth-order-a c curate mesh-adapted differencing 
(|Baker fc van Meteil l2005l ). Interpolation between 
refinement regions was fifth-order-accurate. 

To explore the parameter space of non-precessing spin 
configurations with some mass-ratio dependence, we have 
performed simulations for seven different data sets of 
black holes with unequal masses, as well as one equal- 
mass case. The initial-data parameters for these eight 
data sets are given in Table [1] The spins in these sim- 
ulations were always orthogonal to the orbital plane. 
Our choice of initial tangential momenta was informed 
by a quasi-circular p ost-Newtonian (FN) approximation 
(jPamour et al.ll2000f) . 

Numerically, we can only directly specify the "punc- 
ture" masses mip and m2p of the two holes; to determine 
each hole's physical (horizon) mass m, we first locate 
its apparent horizon (usin g an adapted ver sion of the 
AHFinder Pirect code; s ee Thornbur^ l2003f l . and then 
apply the IChristodoulod (|1970( ) formula: 

— mfj.^ + J^/4m^^^, (1) 

where J is the magnitude of the spin angular momentum 

of the hole, rriirr = ^/AahJIGtt, and Aah is the area of 
the apparent horizon. 

The parameters relevant to our discussion - the mass 
ratio q = TOi/m2 and the dimensionless spin parameters 
ai = oi/toi and 0,2 = a2/m2 - are listed in the first eight 



TABLE 1 

Initial data parameters. Runs are labeled "EQ" for equal 
mass and "ne" for unequal mass. jl and j2 are the spin 
angular momenta of the two holes, either aligned (positive) 
or anti-aligned (negative) with the orbital angular 

MOMENTUM, rripi AND mp2 ARE THE DIRECTLY SPECIFIED PUNCTURE 
MASSES OF THE HOLES. P IS THE INITIAL TRANSVERSE MOMENTUM ON 
EACH HOLE, WHILE L IS THE INITIAL COORDINATE SEPARATION OF THE 
PUNCTURES. 



Run 




J2(M2) 


mpi(M) 


mp2{M) 


P (M) 


L (M) 


NE 


-0.032 


-0.072 


0.374 


0.586 


0.119 


7.05 


NE_+ 


-0.032 


0.072 


0.374 


0.586 


0.119 


7.05 


NEo- 


0.000 


-0.072 


0.374 


0.586 


0.119 


7.05 


NEoo 


0.000 


0.000 


0.382 


0.584 


0.119 


7.05 


NEo+ 


0.000 


0.072 


0.374 


0.586 


0.119 


7.05 


NE+_ 


0.032 


-0.072 


0.374 


0.586 


0.119 


7.05 


NE++ 


0.032 


0.072 


0.374 


0.586 


0.119 


7.05 


EQ+_ 


0.050 


-0.050 


0.480 


0.480 


0.124 


7.00 



rows of Table 121 We have striven to maintain g = 2/3 in 
all of our unequal-mass simulations. 

The grid spacing in the finest refinement region, 
around each black hole, was hf = 3M/160, with the 
exception of our nonspinning unequal-mass case, which 
was on e of a set of runs described previously (| Baker et al.l 
l2006cl ). and for which we used hf = M/40. The ex- 
traction radius was at i? = 45M in every case except 
for the nonspinning case, where it was at i? = 50M. 
For one of our new physical configurations (NE_|_+) we 
also extracted at R — 40M, finding a final kick within 
O.Skms""'^ of that extracted at i? = 45M. Assuming a 
radially dependent error tha t falls off as 1 / R, as fou nd in 
a similar kick computation (Gonzalez et al.' 2007b) , this 
implies that the kick extracted at i? = 45iV/ is within 
8% of what would be computed at infinite radius. Also 
for this physical configuration we ran a higher resolu- 
tion, hf = M/64, to verify that the lower resolution of 
hf = 3M/160 would be sufficient. We found satisfactory 
convergence of the Hamiltonian constraint (Fig[T]) and 
consistency of the radiated momentum (Figl2|) . 



The t hrust dP'^/dt imparted by the radiation was de- 
rived bvlNew man fc Tod! (|1981f) . and is computed from a 
surface integral of the squared time-integral of the radia- 
tive Weyl sc alar ■04 times the unit radial v ector, as given 
explicitly by iCampanelli fc Loustol (|1999D : 




dt 



(2) 



To perform the angular integration in ([2]), we use 
the se cond- order Misne r algo rithm described in iMisneH 
(|2004f) and iFiske eFall (|2005f) . This result is then inte- 
grated numerically to give the total radiated momentum 
P'; to obtain the final velocity of the merged remnant 
black hole, we divide by the final black-hole mass, as 
computed from the difference of the initial ADM mass 
and the total radiated energy. 

3. RESULTS AND DISCUSSION 

In Fig. 121 we present the aggregated recoil kick from 
each of our simulations. The kicks obtained range from 
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Fig. 1. — LI norm of the Hamiltonian constraint for the case 
of unequal masses with prograde spins (NE-|-+) at two different 
resolutions, computed on a domain extending to \xi\ = ASM, which 
includes our wave-extraction surface. The finest level, containing 
nonconvergent points near the puncture, has been excluded. The 
lower resolution, hf = 3Af/160, has been scaled such that for third- 
order convergence it should superpose with the hf = M/6A curve. 




400 



Fig. 2. — The total radiated momentum for the case of unequal 
masses with prograde spins (NE_|_+) at two different resolutions. 
The curves have been shifted slightly in time to line up the peaks, 
for comparison of the merger dynamics. The agreement is better 
than 5% for most of the evolution and better than 2% through the 
merger and final kick {t > 300M). 

~ 60kms~^, in the case where the larger hole's spin is 
aligned, and the smaller anti-aligned, with the orbital an- 
gular momentum, to ^ 190kms~^, when the alignments 
are reversed. The kicks possess a common profile, with 
the bulk of the momentum radiated over ~ 50 M before 
merger. In all but the equal-mass case (EQ., ), we ob- 
serve a sharp monotonic rise in kick over 40 M, followed 
by a substantial "un-kick". That is, around the time of 
merger and ringdown, we often observe a sudden thrust 
in momentum that is directed counter to the momen- 
tum that had accumulated during inspiral. In the EQ^ 

case, this un-kick is absent. We summarize the final kicks 
for each of our configurations in the first eight rows of 
Tabled 



NE_ 
.... NE^ 

- NEo. 
-- EQ^. 

B ■□ ™M 
^ [- NE_^_| 

>^*NE„^ 
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Fig. 3. — Aggregated recoil kicks from all runs listed in Table [T] 
The merger time for each binary matches the peak in its kick pro- 
file; the relative delay in merger times between data sets differing 
in initial spins is consistent with the results of ICampanelli et al.l 
pOO tjc) . AH configurations show a marked "un-kick" after the 
peak, with the exception of the equal-mass case, EQ^ . 

can significantly extend or reduce the time to merger, de- 
pending on whether the spins are aligned or anti-aligned 
with the orbital angular momentum. We observe a sim- 
ilar trend in merger times, as illustrated in the peak of 
the aggregated recoil kicks. This tendency had been also 
been expected based on PN calculations which show that 
the last stable orbit is pushed to s maller radius, imply- 
ing later merger, for aligned spins (jPamou 3 12001 . Note 
that, although resolution can also affect merger time, for 
these short runs we have sufficiently resolved the dynam- 
ics that numerical error in merger time appears negligible 
compared to the effect of spin, as demonstrated in Fig. [51 
Our simulation results, together with data reported by 
other groups, allow us to consider a simple description of 
the total kick for arbitrary mass ratio when the two black 
holes spins are aligned or anti-aligned with the system's 
orbital angular momentum. 

("Herr mann et al.l 



Several re c ent papers ("Herr mann et al.l l2007al : 
i Koppitz et"all 120071 : iChoi et al. 2003) have suggested 
that the kick velocity resulting from comparable-mass 
binary black hole mergers may be approximately 
described in terms of a simple scaling dependence 
consistent with the scaling in the leading-order post- 
Newtonian approximation treatment (Kidder 1995). 
For kicks generated by non-spinning black holes of 
unequal masses, this produces the Fitchett scaling, 
which provides a reasonable approximation to recen t 
numerical simulation results (i Gonzalez et al.l |2007a'). 
Recent papers on spinning black hole mergers suggest 
that the kicks from nearly equal-mass mergers may 
scale linearly, through the quantity A = qai — 02, with 
spins that are either aligned or anti-aligned with the 
orbital an gular momentum (jHerrmann et all l2007at 
lKoppitz"e t al. 2007). For head-on collisions. IChoi et al.l 



It has been noted bv lCampanelli et al.l ()2006cf ) that the 
presence of spins on black holes in an inspiraling binary 



( 200l " lave shown that the PN expressions describe the 
scaling of the spin-asymmetry kick and its relation to the 
kick induced by mass-asymmetry, correctly predicting 
the relative directions of the two kick contributions and 
supporting the idea that the effects of asymmetries in 
spins and masses can be considered independently. 

For inspiraling mergers, with the assumption that the 
radial velocity is small compared to the tangential veloc- 
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ity, the PN prediction for the cases we consider is that the 
instantaneous thrust generated by the spin-asymmetry 
will be aligned with that of the mass-asymmetry, though 
the relative size of these effects may vary as the merger 
proceeds. Unlike the head-on collision case, the direction 
of thrust should vary as the system revolves in inspiral- 
ing mergers, so that we cannot infer that the net effects 
of spin- and mass-asymmetry will produce coUinear con- 
tributions to the overall kick. Motivated by these ob- 
servations we will consider our data with the hypothesis 
that the magnitudes of the kicks induced by spin- and 
mass-asymmetries each scale independently with the PN- 
predicted scaling, but that the directional alignment of 
these two contributions to the kick may differ by some 
angle 9. The total kick may then be of the form 



TABLE 2 

Predicted versus compute d kick speed. Runs label ed 

"S0.##" ARE TAKEN FROM IHeRRMANN ET AL.I (j2007AD . 
WHILE RU NS LABELED "R# " ARE TAKEN FROM 

IKoppitz et AL.I (j2007l 1. 



X -qY + 2{l-q)K COS 6 + , (3) 

where K = k{qai — 02). The parameter Vq gives the 
overall scaling of the kick (note that the factor in brack- 
ets becomes unity for g = 1), while k gives the relative 
scaling of the kick contributions from spin- and mass- 
asymmetries. This expression amounts to a general- 
ization of the post-Newtonian -ins pired kick forrnula dis - 
cussed bv lFavata et al.l (|2004f) and lKoppitz et al.l (|2007[ ). 
which was consistent with coUinearity of the two kick 
contributions, 9 = 0. 

We have tested the simple formula ([3]) with our 
own kick speeds, to gethe r with published kicks from 
IKoppitz et all (|2007D and iHerrmann et all (|2007aD . for 
a total of 17 independent data points. Without more 
precise knowledge of the uncertainties for each measure- 
ment it is not possible to do a true statistical fit. As 
a substitute, however, we have obtained values for the 
three free parameters Vo, 9, and k by minimizing the to- 
tal of a X^-like quantity OC I]i(^pred,i - Wnum,i)^/l'imm,i, 

where Wprcd,! is the predicted, and I'num.i is the measured, 
kick speed for the ith combination of parameters. The 
proportionality constant is chosen so that the minimum 
of "x^" is 14, equal to the number of degrees of freedom. 

With this procedure, our best fit to all simulations cur- 
rently reported gives Vq = 276km s~^, 9 = 0.58 rad, and 
k = 0.85. The minimal regions containing 95% of the 
probability for each parameter are Vq = 267— 294km s~^, 
9 = 0.45 - 0.65 rad, and fc = 0.8 - 0.89. In Table [2] we 
compare the predictions of this formula to our results and 
those of other groups who have explored prograde or ret- 
rograde spins. We see that this simple formula performs 
well, with errors less than ~ 10% in all cases. 

We can estimate the uncertainty in 9, and judge how 
strongly we can rule out a constant = 0, by mak- 
ing the conservative assumption that all the numerical 
kick results have 10% statistical errors. A chi-squared 
analysis then indicates that at one standard deviation 
9 = 0.58 ± 0.8 rad. We also find that, formally, 9 = 
is ruled out at > 5cr and 9 = t:/2 (sum in quadrature) 
is ruled out at > 12cr, but this far from our best values 
the error contours are clearly non-Gaussian. It is pos- 
sible that a more complicated model (e.g., one in which 
9 depends on the mass ratio) is a better representation 
than 0=constant, but no motivation for this exists in the 
current data. 

The success of our fit in describing the existing data 
suggests that this simple expression may, for many astro- 



Run 


q 


di 


02 




1 A^;| 

'^pred Vji^jj 


-{%) 


NE 


0.654 


-0.201 


-0.194 


116.3 


119.5 


2.7 


NE_ + 


0.653 


-0.201 


0.193 


58.5 


58.2 


0.5 


NEo- 


0.645 


0.000 


-0.195 


167.7 


153.1 


8.7 


NEoo 


0.677 


0.000 


0.000 


95.8 


98.6 


2.9 


NEo+ 


0.645 


0.000 


0.194 


76.9 


71.7 


6.8 


NE+_ 


0.655 


0.201 


-0.194 


188.6 


181.9 


3.6 


NE++ 


0.654 


0.201 


0.194 


83.4 


92.4 


10.8 


EQ+_ 


1.001 


0.198 


-0.198 


89.8 


92.6 


3.2 


SO. 05 


1.000 


0.200 


-0.200 


96.0 


93.8 


2.3 


SO.IO 


1.000 


0.400 


-0.400 


190.0 


187.6 


1.2 


S0.15 


1.000 


0.600 


-0.600 


285.0 


281.5 


1.2 


S0.20 


1.000 


0.800 


-0.800 


392.0 


375.3 


4.3 


rO 


1.000 


-0.584 


0.584 


260.0 


274.0 


5.4 


rl 


0.917 


-0.438 


0.584 


220.0 


220.8 


0.3 


r2 


0.872 


-0.292 


0.584 


190.0 


178.1 


6.3 


r3 


0.848 


-0.146 


0.584 


140.0 


141.9 


1.4 


r4 


0.841 


0.000 


0.584 


105.0 


110.4 


5.1 



physical simulations, adequately describe the dependence 
of the kicks on the portion of mass-ratio and spin param- 
eter space that we have studied. However, recent simu- 
lations have suggested that the dominant component of 
the kick may be out of the orbital pl ane deriving from 
spins which lie in the orbital plane (jCampanelli et al] 
llOOT^), a conf iguration outsid e the pa rameter space we 
have studied. iGonzalez et al.l ()2007bl ) have shown that 
such a configuration can produce kicks which may ex- 
ceed 2500 km s^^ directed out of the orbital plane. 
These results appa rently confirm the early predictions of 
iRedmount &: Reei (il989) that out-of-plane kicks would 
be particularly significant. 

Since initial submission of this paper. [Campanelli et al.l 
(|2007a|) have suggested a combined formula, which gen- 
eralizes our formula Q to include out-of-plane kicks 
(u|| in their notation) of the kind discussed above. If 
kick speeds > 1000km s^^ are common in comparable- 
mass mergers of black holes with substantial spin, this 
comes into apparent conflict with the observation that 
essentially all galaxies with bulges appear to have super- 
massive black holes in their cores, sin ce galactic escape 
speed s tend to be < 1000km s^^ fsee iFerrarese fc FordI 
I2OO5I for a review of supermassive black holes and their 
correlation with galactic properties). It seems unlikely 
that most supermassive black holes have low enough 
spins to guarantee small kicks, given evidence such 
as broad Fe Kg line s from a number o f black holes 
llwasawa et al.' '1996'; 'Fabian et al.' '2002": 'Miller et alj 
j|op2- Reynolds & Nowak 2003; Reynolds et al. 2003 
iBrenneman fc Revnoldsl 120061) as well as overall argu- 
ments from the inferred hig h average radiation efficiency 
of supermassive black holes (|Sohanlll982l : lYu fc Tremaind 
l2002f ). iBogdanovic et al.l ( 20071 ) suggest that torques 
from gas-rich mergers tend to align the spins of the holes, 
and thus lead to small kicks, but no known preferred 
alignment exists for gas-poor mergers. More exploration 
of the parameter space of spin kicks is clearly necessary. 
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